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Abstract 



We have observed multi-photon resonances in a system with a spin 3/2 ir- 
radiated simultaneously by a multiple pulse radiofrequency sequence and a 
low frequency field swept in the range 0-80 kHz. The used excitation scheme 
allowed us to measure the effective field of the radiofrequency sequence. A 
peculiarity of this scheme is that the intensity of the resonance lines decreases 
slowly with the mode number. The theoretical description of the effect is pre- 
sented using both the rotating frame approximation and the Floquet theory. 
Both approaches give indentical results at the calculation of the resonance 
frequencies, transition probabilities and shifts of resonance frequency. The 
calculated magnetization vs. the frequency of the low-frequency field agrees 
well with the obtained experimental data. The multi-frequency spectra give 
a way for studying slow atomic motion in solids. 
Introduction 

One of the most effective and promising high-resolution nuclear magnetic resonance 
(NMR) and nuclear quadrupole resonance (NQR) techniques for the study of solids is a 
multiple-pulse radiofrequency (RF) action [1,2]. The multiple-pulse methods allow one to 
remove dipolar broadening from a resonance line in solids thus increasing by several orders 
the sensitivity of the NMR and NQR spectroscopy in the study of weak interaction. These 
methods are very effective in the study of the spin lattice relaxation processes due to a slow 
atomic motion. Usually the theoretical description of multiple-pulse experiments both in 
NMR [3] and NQR [4] is based on the construction of the effective time-independent Hamil- 
tonian by using the conditions for periodicity and cyclicity of the pulsed action . Then the 
dynamics of a spin system subjected by pulsed RF fields is presented in an equivalent form 
as the motion of nuclear spins in a constant effective field H e [5] . The magnitude and direc- 
tion of this effective field are determined by parameters of the multiple-pulse sequence. An 
experimental measurement of the value of the effective field is important for the confirmation 
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of this theoretical model. 

It is reasonable to suggest that an additional field with an angular frequency Q close 
to uj e = should cause resonance absorption of energy (7 is the gyromagnetic ratio of 
nuclei). Spin-echo signals observed between RF pulse sequence would allow us to determine 
H e as well as to obtain the information on slow atomic motion that is not available from 
the traditional metrods 

With this in mind, we have studied experimentally resonance transitions in the nuclear 
spin system subjected by a simultaneous action of a multiple-pulse RF sequence and an 
additional low frequency (LF) field with an angular frequency fl. The results of our exper- 
iments described in the next section have shown that resonance transitions were observed 
not only at the frequency close to Q = u e , but also at frequencies close to Q n given by the 
expression: 

tt n = \u e ± 2im/t c \ , n = ±l,±2... (1) 

where t c is the period of the multiple pulse RF sequence. 

Multiple resonance modes of higher orders have been detected by microwave spectroscopy 
[6], molecular beam technique [7], optical pumping [8] and observed previously in NMR 
experiments [9-13]. However the amplitude of these resonances decreased abruptly with 
the mode order of the resonance. As distinct from this, the amplitude of the resonances 
observed in our experiments decreased slowly and the resonances of higher orders were well 
observable. 

Because the nuclear spin system possesses a set of the resonance frequencies, the relax- 
ation measurements performed on one resonance frequency Q n can give the information on 
oscillations of atoms on all the frequencies from the spectrum determined by (1). It moves 
us to comprehensive experimental and theoretical study of this system. 

The theoretical treatment of NMR phenomena is usually based on three approaches: i) 
a semi - classical mathematical approach [14]; ii) a second quantization method [15]; and 
iii) the Floquet theory [17]. 
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The semi - classical mathematical approach [14], where the field is considered as a classical 
system and the atomic system as a quantum one, has allowed one to explain of series of 
experimentally observed phenomena. This approach is quite natural if to take into account, 
that the average number of photons in a mode of the periodic field is extremely great. 
The main method used in the framework of the semi - classical approach is the so-called 
"rotating frame approximation", keeping exactly just the terms that are resonant. The 
remaining non-resonant terms are considered as a perturbation. 

Intrinsic inconsistency of the semi-classical approach is obviated in the framework of the 
secondary quantization method [15,16]. Treating the RF field as photons, the evolution of 
the united system "atom +field" (so called "dressed" atom) is described by the Hamilto- 
nian which is independent of time, and its investigation turns out simpler than solving the 
Schrodinger equation with the time-dependent Hamiltonian. With the time-independent 
Hamiltonian, one can define energy levels of the physical system. Each of these levels cor- 
responds to a stationary state of the system: atom dressed by photons. The dependence of 
the energy of the united system on the field frequency allows one to construct the diagram 
of levels. The coupling between RF photons and an atom perturbs the energy levels and 
allows one to interpret all resonances as "crossings" and "anticrossings" of the conforming 
levels [18]. 

The Floquet theory [17] is a power method widely used in NMR spectroscopy for solving 
time-dependent problems [10,19]. On the one hand, this approach uses the advantages of 
the secondary quantization method resulting in the time- independent Hamiltonian. On the 
other hand, it simplifies calculations as the terms related to the system of free RF photons 
are not included in the Hamiltonian . The last is justified because an RF field applied in the 
NMR method contains abundant photons and changes in the state of the photon system may 
be neglected. In the framework of the Floquet theory, the interaction of the spin system with 
RF field is considered as completely quantum-mechanical. Therefore, the Floquet theory 
can be considered to be a bridge between the semi-classical and the quantum methods. 

To calculate the transition probabilities and the shift of the center of the resonance line, 
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two theoretical methods are used: the semi-classical approach and the Floquet theory. We 
extend the theory of the multiphoton resonances [14,15,17] to the case of the pure multiple- 
pulse NQR. The results given by both theoretical approaches are compared with each other 
and with the experiment. 
Experiments 

The experiments were performed using an automated multiple-pulse NQR spectrome- 
ter. Mutually perpendicular continuous LF and pulse RF magnetic fields were generated 
by crossed coils. Resonances were observed in the effective field of a multiple-pulse RF 
sequence (7r/2) — (t c /2 — <p — t c /2) N , where ip denotes the pulse that rotates the nuclear 
magnetization about direction of RF field in the rotating frame by an angle ip, and N is 
the number of pulses in the sequence. This sequence consisted of N = 256 pulses and spin 
locking signal was sampled in the interval between them. The period of the pulse sequence 
was t c = 100ms. 

The NQR of 35 Cl nuclei was observed in polycrystalline KCIO3 at 77 K. Without action 
of the LF field, we recorded the sequence of echo signals in the interval between pulses at 
28.9539 MHz with practically constant amplitude (height of the spin echo signals) during 
the observation time (Fig. 1). This amplitude corresponds to a quasi-equilibrium state of 
the spin system with magnetization Mq. 

LF field with the amplitude 2.5 G, nonsynchronized with RF pulses, was swept in the 
range 1 -j- 80 KHz. The resonance reduction of the magnetization amplitude was observed 
(Fig.l) at several different frequencies Q . The transient process with the spin-spin relaxation 
time T 2 was followed by establishment of a new quasi-equilibrium state of the spin system 
with the reduced amplitude of the magnetization. Our measurements give the value T 2 =455 
ms. 

Fig. 2 shows the dependence of the relative magnetization M/M on the frequency Q of 
the LF field for the ip = f RF pulse (corresponding to the pulse duration t w — 15 /is ). The 
effective frequency of the RF field is determined by u e = j- . For example, u e = 2.5 kHz 



5 



for ip = f . 

As follows from Fig. 2, the amplitude of the resonances decreases slowly with increasing 
the mode number n. 
Theory 

1. Semi-classical approximation: rotating frame approximation 

Let us consider a system of nuclear spins 1 = 3/2 placed in an inhomogenous electric 
field gradiet (EFG) and subjected to a joint action of two time-dependent magnetic fields: a 
multiple-pulse RF sequence with the angular frequency to equaled to the resonance frequency 
u and a continuous low frequency (LF) field with the angular frequency Q. 

In the rotating frame the equation for the density matrix of the system takes the form: 

dp (t) 



dt 

where the Hamiltonian of the system is 



= [H(t),p(t)], (2) 



H (t) = 2lu 2 S z cos ttt + S x cp^2S(t- kt c - t c /2) , (3) 

k=l 

S x and S z are X- and Z-components of the effective spin operator [4]; the pulse angle 
(p = ^Hity, and < <p < 2n; cv 2 = lH 2 ; Hi and H 2 are the amplitudes of the RF pulse and 
LF fields, respectively. The initial phase of LF field is chosen equal zero. 
To solve Eq.(2), we apply the unitary transformation 

p {t) = e ~ ms * t P+ (t) p (t) P it) e ins *\ (4) 

where the unitary operator P (t) is given by the solution of the equation 

= S xV jr5(t- kt c - t c /2) P(t)-P (t) H e (5) 
fc=i 

with the initial condition 

P(0) = 1 (6) 
and the effective time-independent Hamiltonian 
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He = UJ e S x . (7) 

After performing the unitary transformation (4) we obtain the following equation for the 
transformed density matrix p (t) 

dp {t) 
dt 

with 



[n tr (t),p(t)} (8) 



H tr (t) = (u e — Q)S X + S z (t) , (9) 

where S z (t) = P + (t) S Z P (t) is the periodic function of time and can be expanded into the 
Fourier series 

oo oo 

S z (t) = S + Ke~ M + S- 6 »e mWc * (10) 



n=—oo 



where S± — S z ± iS y , b n = ^2 J + 2™ n 2 are ^ e F° ur i er coefficients, and u c = j 1 . 

The resonance terms in the solution of Eq. (8) are determined by the lowest values of 
the differences £ = \u e — Q ± nu c \. At £ = we obtain the expression (1) for resonance 
conditions. The rest non-resonant terms determine a frequency shift A for the position of 
the resonance dips shown in Figs 1 and 2. Let us estimate this shift as a function of the 
frequency Q using the average Hamiltonian theory [24,21]. We will consider three various 
cases depending on the value of the frequency Q. 

The first case, Q = 0, corresponds to the appearance of an additional constant magnetic 
field along Z-axes resulting in the difference between u and ujq. 

To determine the influence of all the non-resonant terms, we will perform the unitary 
transformation under Eq.(8), 

f /(t)=e ivS '-hp(t)e- ivS '-h, (11) 

which leads to 
dp' (t) 



dt 



[n'(t), P '(t)} 



00 



S+ b n e- l(27Tn+ip) - + S- b n e (27Tn+ip) -,p'(t) 



(12) 



Let us replace H' (t) by the average Hamiltonian up to second order of the expansion in 

u 2 /uj e : 

°° b 2 

n=— 00 

Note that in the case of non-zero initial phase of the LF field H includes \b n \ 2 . Therefore 
the initial phase of the LF field does not influence the result. Performing the summation in 
Eq.(13) we obtain 
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Ai = — cot -ip with ip 7^ — 2-kti. (14) 



In the limit that t c — > 0, a one long pulse, we obtain from Eq. (14) a simple expression 

2 

for the shift A 1 = which is similar to the Bloch-Siegert results [22]. However, in contrast 
to the Bloch-Siegert shift which is produced by the RF field and has the value of the order 
of (jj^ ~ 1CT 6 , the obtained shift is sufficiently higher because it is determined by the 
ratio of two weak fields with the amplitudes H 2 and H e . For example, for the angle ip = f , 
the shift Ai is of the order of 10~ 2 . The shift Ai is caused by a deviation of the frequency 
of the applied RF field from the Larmor frequency and should be taken into account when 
the results of multiple-pulse experiments are analyzed. 

The second case: Q = u) e . The average Hamiltonian H is 

00 b 2 

H av =2b S x + V (15) 
A — ' , lnuj c 

n=—oo,riyL0 

After the summation in Eq.(15) we obtain the shift due to common action of the non- 
resonant terms 



2^ „.„2 1 .. <P 



If 3 



sin 2^ _ 8 ^ + Sin ^ 



(16) 



Thus, the resonance frequency differs from uj e . A shift appears also when we consider 
Q = u e + nu c . 

In order to determine correctly the resonance conditions, let us consider the third case: 
VL 7^ u e + nu c . Using the transformation 



P" (t) = e 



tc 



p(t)e 



-i(<p-$)S x 



(17) 



where $ = Qt c , we obtain 
dp" (t) 



dt 



n=—oo n=— oo 

Let us exclude from the sums in (18) the terms with n = — : 
dp" (t) 



(18) 



dt 



(-l)^^sin|^ + S+ 

_j_ tp — Q 
1 1 Ztt 



^2 bnC 



—i(ip—$+2-Kn)- 



+s- E 



b n e 



i((p-<J>+27m)-^ n ii 



n=~oo,n^-^— 



(19) 



The average Hamiltonian up to second order of the expansion in u; 2 /cj 6 
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(if + 2nn) (</?-$ + 2?™; 



[S-,S + ]. 



(20) 



Because in the considered case f — $ ^ —2nn, the first term in (20) disappears and the 
summation over all n gives the following equation for the frequency shift: 



A 3 = u 2 2 t c | 1 



2 sin | sin ^ 
$sin^ 



(21) 



Emphasize that $ is a function of A 3 and Eq. (21) has infinite number of roots. Reso- 
nance transitions are realized at fl = Cl n , where 



VL n = u e + A 3n + nuj c , n = 0, ±1, ±2, ... 



(22) 



Resonance frequencies calculated according (22) using numerical solutions of (21) for the 
shift are presented in Table 1 and 2 along with the experimental data. The shift is caused 
by off-resonant component of the LF field. 

The time-average transition probability, P can be determined using the method devel- 
oped in [23] : 



ME 



2^ (u e -n-ncu c + A 3 ) 2 + bl' 



(23) 



Since « 1 (fc> r example, = 6. 5 x 10~ 2 for tp = |), the transition probabil- 

ities become significant only in the vinicity of the resonance frequency determined by the 
condition (22). The expression for the nuclear magnetization M , the quantity needed to 
compare with experiment, can be found using Eq. (23): 

Mo ^ (u e -n-nu c + A 3 f + bl 

Expressions (23) and (24), obtained in the framework of the semi-classical approach, deter- 
mine a multi-line resonance spectrum of the two-energy level system with spin 1/2 irradiated 
by multiple pulse RF and LF fields. As follows from (24) the amplitudes of the resonances 
decreases slightly with the number n. There is disparity between a two-level energy model 
of the spin system and the multi-line resonance spectrum. In order to make a physical 
interpretation clear, the energy spectrum and transition probabilities will be calculated in 
the next section by using the Floquet theory [17]. 

2. Floquet theory 

The Floquet theory provides a quite natural description of the multi-line resonance spec- 
trum of the two-energy level spin system. Let us introduce the time-dependent operator 

U (Mo) = F (t) F + (t ) , (25) 

where the unitary operator F (t) obeys to the evolution equation 

H tr (t) F (t) (26) 



dt 

with the periodic time- dependent Hamiltonian (9). The general form of the solution of 
Eq.(26) is given by the Floquet theorem [17] is 

F (t) = * it) e- mF \ (27) 

where \& (t) is a periodic function of time with the frequency u c , and 7i F is the time- 
independent effective Hamiltonian Using the periodic character of the Hamiltonian H (t) 
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and function ^ (t) , the time-dependent Eq. (26) can be rewritten as an infinite set of the 
coupled equations [17] 

J2 [n n a ' k + nu c 5 nk 5 a ,) = hpFZp, (28) 

fik 

where hp are diagonal elements of the effective Hamiltonian H F , 7i"~ fe and are the 
Fourier components of the matrix elements of H* (t) and F (t) : 

F a p{t) = Y,F^ mWct z- lhp , (29) 

n 

Kp(t) = J2 n «P ei ™ ct - ( 3 °) 

n 

Here a = ±1/2, j3 = =Fl/2 characterize various spin states. Using the orthogonal basis \an) 
, where n denotes the Fourier component, the term in the brackets of the Eq. (28) can be 
presented as the Floquet Hamiltonian [17] 

(cm| H F \0k) = HI/ + ruv c 5 nk 5 a p, (31) 

Using the operator basis e an ^ m [24] with matrix elements (/ip\ e an ^ m \qq) = 5 aix 5p v 5 pn 5 mq 
the Floquet Hamiltonian can be expressed as 

7~(-F = ^ ] ( £ q + niO c ) e an ^ an + ^ ] (bji-n e an,f3m + C-m e /3n,Qtnj (32) 
a,n a^P,n,m 

where e a = a (cu e — Q) , b°^_ n = b m ^ n are the Fourier coefficients determined in (10). The 
first term in (32) is diagonal, while the second one is off-diagonal. 

As follows from (32), nonzero off-diagonal matrix elements between \an) and \/3m) for 
various n and m appear in first order of Therefore all resonance transitions between 
different states of the Hamiltonian (32) can be excited in the straight way without involv- 
ing intermediate states as it is used in the scheme of multiphoton excitation proposed in 
[10,19,25] . 

The intensity of a resonance transition can be estimated using the perturbation method 
[7] based on the approximation of the Floquet Hamiltonian Hp (32) by the 2x2 matrix 

(a ^(3) 
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£ a -\- TTLU) C -\- c am b m — n -\- d n _ m 



^m-n + d n -m £ f3 + ™c + Cp n 

with matrix elements corrected by taking into account the rest part of the Hamiltonian: the 
diagonal elements - by introducing the additional terms c an which are up to second order to 
within (^j < 1 is 

E(an\ Tip \fJ>k) (/ik\ TIf \om) 
E —Ei, ' 



Can 



and the off-diagonal ones - by adding the terms d°^_ m of the same order 

jq/3 =2 y (Pn\H F \fJ.k) (fJ.k\H F \arn) ■ , . 

In (34) and (35) -Eg n are eigenvalues of the Floquet Hamiltonian H F . Diagonalizing the 
matrix (33) gives the eigenvalues of the Hamiltonian H' F : 

2E%? ) = [e a + e f) + (m + n)u c ]± 

± { [(e a + muj c + c am ) - (ep + nu c + cp n )] + b 

• (36) 

Substituting the matrix elements (31) into (34) we obtain the expressions for the energy 
level shifts 

°° b 2 

Can = -Cfin = (y, _ $ + 2tf (jfc - 7l)) ' (37) 

The summation in Eq. (37) gives the same equation for the frequency shift as above obtained 
using the rotating frame approximation (see Eq.(20)): 



LVot r ( . 2 sin | sin * 
$ sin 



A, = c„„-c,„=^|l- .~'~' |. (38) 
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As follows from Eq. (35), the terms d°^_ m are proportional to {^f^ and can be neglected 
in Eq. (36). Therefore the energy levels are given by 

E%f> = \ {(m + n) u c ± [{u e - Q - (m - n) u c + A 3 ) 2 + b^ m ] 1/2 } . (39) 
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The normalized energy — is presented in Fig. 3 as a function of — for m — n = 0. 




The average probability [17] of the resonance transition is 




2 ("e " « " (m - n) u c + A 3 ) 2 + fc^_ n ' 



(40) 



The probability calculated according to Eq. (40) is presented in Fig. 4 at tp 
of x = — for the terms in the sum with m — n = 0, ±1, ±2. 



I as a function 



Using (40) the expression for the nuclear magnetization M, is the quantity needed to 
compare with experiment, can be found: 



Eqs. (40) and (41) coincide (by changing the summation index: k = m — n) with the 
equations obtained for the transition probability and the magnetization in the previous 
section (see Eqs. (23) and (24)). The magnetization calculated according to Eq. (41) is 
presented in Fig. 5. 

Results and Discussion 

As follows clearly from the comparison of expressions (23), (24) , (40), and (41), the 
calculations using both the semi-classical and Floquet methods give the identical results 
for the transition probability and nuclear magnetization. The results for the resonance 
conditions and for the shift A 3 of the resonant frequency obtained within the framework 
of these methods are also identical. Moreover, the semi-classical approach allows one to 
obtain automatically the shift Ai under off-resonance conditions uq ^ uo {Vt = 0) from 
RF excitation. However the interpretation of the multi-frequency resonance transitions for 
two-level energy system in the framework of the semi-classical approach involves difficulties. 

The Floquet method gives clear physical description of observed phenomena. The time- 
independent Hamiltonian obtained using this method leads to the multi-level energy spec- 
trum. Each of these levels corresponds to a stationary state of the system: spin dressed by 
photons. The coupling between photons and spin allows one to interpret the resonances as 
" anticrossing" of the levels [17,10]. 




(41) 



13 



Matching of the computed results for the nuclear magnetization with the experimental 
data displays series of essential differences. One of these differences consists in the fact that 
the observed NMR signal does not damp up to zero (Figs. 1 and 2), while the calculations 
predict zero magnetization at the resonant frequency (Eq.(41) and Fig. 5). A possible reason 
of the observed finite value of is an inhomogeneity of the RF field in the bulk of the 
sample resulting in a distribution of angles (p over the sample. To compare correctly the 
theoretical result for the nuclear magnetization with the experiment, we have to take into 
account the distribution of <p over the real resonance line with nonzero line width averaging 
expression (41) over the line: 

® = s; I"** <">'<">■ (42) 

Let us consider the Lorentzian distribution of frequencies over a line having width c . It 
leads to the following distribution of angles (p 

g(x) = -^— 2 , (43) 

where x = ip — ip and ip = (<f) v is an angle averaged over the sample. After the integration 
of (42) with the substitution of (43) we obtain 



M\ =1 _£ l + a^_ (M) 



Mi 



where 



1 + — n , (45) 



(cu 2 t c ) 2 sin 2 f V V?o 
The average magnetization, calculated according Eq. (44) with c = 0.018, is presented 



in Fig. 6 at ip = | as a function of x 



along with the experimental results. Using the 
relationship c = ((<p — ipo) /<fo) the dispersion of angles can be estimated as ((</? — ipo)) = 
2.8274 x 10 _2 rad or 1.62°. One can see the excellent agreement of the theoretical and 
experimental results. 
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The resonance transitions corresponding to negative frequencies Q (see Fig. 5) can be 
interpreted as ones caused by the rotating in the opposite direction component of the linearly 
polarized LF field . 

The multi-frequency resonances were observed in the experiments with other excitation 
schemes [10,19,25] and explained by multiphoton transitions between different states. Be- 
cause the multiphoton transitions involves intermediate states and the transition probability 
decreases with number of photons [10,19,25], the intensity of the resonance signal dimin- 
ishes quickly with number of the mode. A peculiarity of the excitation scheme proposed in 
the present paper is that, along with multiphoton transitions, there are direct transitions 
between dressed spin states with (n — m) > 1 . Therefore, the transition probability is the 
same for all the differences n — m. As a result, the intensity of the resonance signal is slowly 
decreased with number of the mode (Figs. 1 and 2). 

Conclusions 

We have studied both experimentally and theoretically the dynamics of a spin system 
with the spin 3/2 under simultaneous influence of DC, multiple-pulse RF and LF fields. An 
important peculiarity of the applied excitation scheme is the possibility of the measurement 
of the effective field H e and observation of odd and even resonances with large mode numbers. 
It was shown that the intensity of the resonance lines decreases slowly. 

Both used theoretical approaches, the semi-classical and the Floquet methods, give the 
identical results at the calculation of the resonant frequencies, magnetization amplitudes and 
shifts of resonance frequency. These results are in a good agreement with the experimental 
data for different angles and mode numbers equaled to n — 0, ±1,±2 . Some difference 
between the theoretical and experimental values of the signal intensity can be explained by 
the inhomogeneity of the RF field over the sample leading to the dispersion of the angle ip . 

The proposed technique can be incorporated for the study of the slow atomic motion in 
systems involving spin-3/2 nuclei. 
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Figure captions: 

Fig. 1 Time dependence of relative magnetization in the absence of the LF field 
(open circle) and in the present of the LF field (solid circle) at ip = | in hexamethylbenzene 
C 6 (CH 3 ) e . 

Fig. 2 Dependence of relative magnetization on the LF field frequency ft at (p — | 
in hexamethylbenzene Cq(CH 3 )q. 

Fig. 3 The normalized energy as a function of x = ^- for m — n = 0. 

Fig. 4 The probabilities as functions of x — — for ip = | and m — n = 0, ±1, ±2. 

Fig. 5 Normalized magnetization as function of x = \fl/u e \ calculated according 
Eq.(41) for <p = \ and m - n = 0, ±1, ±2. 

Fig. 6 Dependence of averaged magnetization of x = |il/o; e | at y?o = | in hexam- 

ethylbenzene Cq(CH 3 )q: solid curve - theory (44), open circle -experiment. 
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